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Abstract 

Problems involving approximation from scattered data where data 
is arranged quasi-uniformly have been treated by RBF methods for 
decades. Treating data with spatially varying density has not been in- 
vestigated with the same intensity, and is far less well understood. In 
this article we consider the stability of surface spline interpolation (a 
popular type of RBF interpolation) for data with nonuniform arrange- 
ments. Using techniques similar to [10) . which discusses interpolation 
from quasi-uniform data on manifolds, we show that surface spline in- 
terpolation on ]R d is stable, but in a stronger, local sense. We also 
obtain pointwise estimates showing that the Lagrange function decays 
very rapidly, and at a rate determined by the local spacing of datasites. 
These results in conjunction with a Lebesgue lemma, show that sur- 
face spline interpolation enjoys the same rates of convergence as those 
of the local approximation schemes recently developed by DeVore and 
Ron in [4] . 

1 Introduction 

Since their introduction as tools for approximation and interpolation in the 
1970s in works of Duchon [6j [T] and Meinguet [H], surface splines (known 
also as polyharmonic splines) have become one of the most prominent exam- 
ples of radial basis functions (RBFs). The hallmark of the RBF methodol- 
ogy is to compute with linear combinations of scattered translates of a basic 



*2000 Mathematics Subject Classification: 41A05, 41A25, 46E35, 65D05 
' Key words/interpolation, surface splines, Lebesgue constant, radial basis function 
department of Mathematics, Texas A&M University College Station, TX 77843, USA. 
Research supported by an NSF Postdoctoral Research Fellowship. 



1 



function (the RBF): 



s = 



5es 



where the linear combination is taken over a finite set of centers E (al- 
though often one supplements a polynomial of low degree to such a linear 
combination). By choosing the RBF correctly, one can effectively treat 
high dimensional data by, e.g., interpolation, least squares minimization, or 
Tikhonov regularization to name a few popular schemes. The surface splines 
are defined by 



In this article, we investigate error and stability of surface spline interpola- 
tion, two subjects which still harbor many open questions. 

Of principal interest are questions concerning norms of the operators 
underlying the approximation scheme, especially when those operators are 
projectors (as is the case for interpolation and L2 minimization), because 
such norms measure the stability of the approximation process. When the 
approximation scheme is interpolation - the scheme investigated in this arti- 
cle - the Loo operator norm is called the Lebesgue constant. Finding bounds 
for such norms based on the computational cost (measured as the dimen- 
sion of the space of approximants) is a pervasive problem in approximation 
theory, despite the fact that there are few examples where such constants 
can be satisfactorily bounded. Univariate spline theory provides an impor- 
tant exception, with operator norms being uniformly bounded; results along 
these lines hold for interpolation, L2 minimization and many other schemes 
- with a bound independent of the complexity of the problem (i.e., regardless 
of the number of knots), cf. [H El I15| . For polynomial interpolation, inter- 
polation by trigonometric polynomials on the circle, or spherical harmonics 
on d-dimensional spheres, Lebesgue constants grow with the dimension of 
the space of interpolants, indicating that interpolation in these settings is 
unstable. We develop and investigate a "penalized" Lebesgue constant, a 
quantity that measures stability of the interpolation process under highly 
nonuniform distributions of centers. 

The results we present in this article show that under some simple as- 
sumptions, surface spline interpolation of compactly supported functions is 
stable, but, more than this, it is stable in a very strong sense, one that takes 
into account the local distribution of centers. 




when d is odd; 
when d is even. 
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In the course of our analysis, we also show that the Lagrange (or fun- 
damental) functions associated with this kind of interpolation are very well 
localized: they decay nearly exponentially in terms of the spatial variable, 
and this decay becomes dramatically more rapid in regions of high density - 
the Lagrange function also scales linearly with the local spacing of the cen- 
ters. This is a significant accomplishment in its own right, as it demonstrates 
that the space of RBFs has a basis that is fairly local. The importance of 
this stems from the fact that RBFs are globally supported and have influ- 
ence far from their center. Even considering RBFs with compact support 
does little to diminish this fact; because one does not dilate the RBF as the 
centers become dense, the RBF captures more and more centers in its sup- 
port as the density increases. The Lagrange functions do, however, adapt 
to changing density of the centers. 

Another result we present concerns the error incurred by interpolation. 
Recently, DeVore and Ron [4] have developed error estimates for surface 
spline approximation that are both precise (with convergence in L p cor- 
rectly dictated by the L p smoothness of the function) and local (with point- 
wise convergence controlled by the local density). In the larger context of 
approximation theory, precision is not so remarkable, but for kernel based 
approximation, precise results have been notoriously difficult to prove. With 
few exceptions, for RBF approximation such results are known only when 
error is measured in L2 , and only then for target functions residing in certain 
reproducing kernel Hilbert spaces. The estimates in [3] are local in the sense 
that the error is also controlled by the local density of the centers. Such re- 
sults, which have long been known for univariate spline approximation, allow 
faster convergence in regions of high density. 

The key arguments used in our results have been developed gradually, 
over a period of more than 30 years. At the heart of our approach is a dis- 
crete Gronwall's inequality employed by Descloux [3f and Douglas, Dupont 
and Wahlbin [5]. A similar inequality was used by Matveev in [13] . to obtain 
results about interpolation by Z) m splines, which are interpolating functions 
satisfying certain variational problems. As in the RBF case, such inter- 
polants can be expressed as a linear combination of kernels, but the kernels 
depend strongly on the geometry of the underlying domain O. Because of 
this, they are not generally well suited for direct computation. But when 
Q = M. d , this corresponds to the problem of surface spline interpolation on 
M. d , and provides an ideal starting point for our setting. Indeed, Matveev's 
work plays an important role in Johnson's |12| resolution of the L p saturation 
order for surface spline interpolation on bounded domains when 1 < p < 2. 

Very recently, this approach has been used to bound Lebesgue constants 



3 



for interpolation on Riemannian manifolds by Narcowich, Ward, and the 
author [TO]. In fact, it is shown in [TO] that when the manifold is M. d and 
centers are assumed to be quasi-uniform, the boundedness of the Lebesgue 
constant for surface spline and 'Sobolev spline' interpolation follows rather 
directly from the results of Matveev |13j . Whereas the challenge in the 
manifold setting is to construct kernels suitable for the type of argument 
developed for D m splines, and then to formulate such an argument to work 
in different geometries, the crucial difference in this article is to find such an 
argument that treats centers with spatially varying density (as opposed to 
quasi- uniform centers whose density can be measured by a single, spatially 
fixed parameter). 

Organization 

In the following section, we introduce the penalized Lebesgue constants and 
describe their role in the stability of the interpolation problem we consider. 
We also discuss the Lebesgue lemma, one of the main motivations for this 
investigation. In Section 3, we discuss the main technical lemmas needed to 
provide sharp bounds for the Lagrange functions. Section 4 provides Sobolev 
and pointwise decay estimates of the Lagrange functions. The main results 
- the bound on the Lebesgue constant and the use of the Lebesgue lemma 
to obtain optimal local approximation rates for surface spline interpolation 
are given in Section 5. 

Background 

In this paper, we treat the interpolation of functions on an open, bounded 
domain C M. d . We denote the diameter of Q by r\ := sup^g^ dist(x, y). 
The functions we interpolate are compactly supported in Q, with 

ro := dist(supp(/), dQ) > 0. 

The set of nodes for interpolation coincide with the set of centers H C !! 
(introduced previously), and this is assumed to be finite. Moreover, the 
subset Ef = Hnsupp(/) is assumed to be nonempty. A consequence of this 
assumption is that the notorious 'boundary effects' are avoided^. 

The space of polynomials of degree less than or equal to n is denoted by 
II n . The set of centers H is said to be unisolvent with respect to LT n if the 
only p £ II n that vanishes on S is the zero polynomial. 

lr The presence of a boundary poses an obstacle in nearly all computational problems; 
approximation with RBFs is no exception to this, although the negative effects of a bound- 
ary and how to overcome these is gradually being understood, cf. [9]. 
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We are now ready to define the surface spline interpolant. Given a 
continuous function / 6 C(R a! ), a set of centers S unisolvent with respect 
to n m _i, and a surface spline cj) = 0m with m > cZ/2, the surface spline 
interpolant to / is the function /=/ that satisfies /=;/(£) = /(£) for all 
(GH, and that resides in the set 

5™ :- ^^(■-5)+p(')|pen m _ 1 , ^A c g(£) = o, v g en m _! 

(1.1) 

Existence and uniqueness of the surface spline interpolant is a consequence 
of the conditional positive definiteness of the function <p m , which, in turn, 
is a consequence of the positivity of its generalized Fourier transform. See 
[16\ Chapter 8] for the relevant background. 

There is a dual, variational characterization of surface spline interpola- 
tion which states that the interpolant In/ is the minimizer of the L2 Sobolev 
seminorm of order m over all interpolants to / at the set H (cf. |16l Chapter 
10]). In other words, \Isf\w^(R) — l u lw 2 m (R d ) f° r an u f° r which u\s = /|s- 
The Sobolev seminorm is defined as follows. 

Definition 1.1. For an open set T C R rf , we define the m th Sobolev semi- 
norm: 



l/k™(T) : = [ E 7ff / l^/(*)l a d^ • (1-2) 

\\/3\=m ' T / 



The space of locally integrable functions having finite m th seminorm is 
the homogeneous Sobolev space, denoted W^(T). When T = M. d , we sup- 
press the dependence on T, writing simply W™, and express the seminorm 
as \f\ m . In the literature this is also frequently referred to as the Beppo-Levi 
space BL m . When m > d/2, the homogeneous Sobolev space is embedded 
in the space of continuous functions, and is, therefore, a reproducing kernel 
semi-Hilbert space. On occasion we also use the (inhomogeneous) Sobolev 
norm \\f\\ w ^ r) := QXa l/l^ (x) ) 1/2 - 

Throughout this article we also consider Sobolev seminorms when T is a 
ball, the complement of a ball, or an annulus. For R > we denote the open 
ball centered at x having radius R by B(x,R). We denote its complement 
by B c (x,R) = R d \B(x,R). For R,w > 0, we define the annulus centered at 
x having outer radius R and width w by 

A(x, R, w) := B(x, R)\B(x, R - w) = B C (x, R - w)\B Z (x, R) (1.3) 
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where equality is modulo sets of measure zero. When x = 0, we will suppress 
the first argument in the above notations. 



Penalized Lebesgue Constants and Surface Spline 
Interpolation 



A motivating goal is to provide estimates of the form 



<C a 



f 



for / G C(fi) 



(2.1) 



where p : O — > (0, oo) is a positive function measuring, roughly, the density 
of the set of centers S: p(x) indicates the spacing of centers near to x. We are 
mostly interested in situations where the density of H is permitted to vary 
considerably. This is done without the usual assumption of quasi-uniformity, 



h := maxdist(x,H) < coq, where q :- 



mindist(£,~ \{£}) 

£6S 



an assumption in effect for most interpolation results; in fact, this is most in- 
teresting when the centers are deliberately in violation of the quasi-uniformity 
assumption. By allowing points to coalesce in such a way that the density 
of points varies spatially, we find that the notion of the global fill distance h 
as an approximation parameter is simply too crude. At some locations we 
will have spacing of points much smaller than h. 

On one hand, an estimate of the form f)2. 1 1) shows the stability of the 
interpolation process: setting a = shows that the interpolant is never 
much greater than the original function, while choosing a > gives a refined 
notion of stability, putting greater emphasis on regions where p is small. 

On the other hand, this leads to a Lebesgue lemma that respects the 
local distribution of centers, one where we can interpolate the error from 
any other surface spline approximant: 



f-Inf 




P° 


oo 



f 



w 



p° 



<d + c 



i inf 



f 



p° 



Here S& is the finite dimensional space of surface spline approximants, de- 
fined in (11. lj) . This shows interpolation is "near best" in the sense of local 
approximation. I.e., that the best local approximation rates are attainable 
for interpolation as well. DeVore and Ron [4] have recently investigated 
local surface spline approximation, providing a scheme / i— )■ st g that de- 
livers precise local approximation results of the form |/(x) — Sf^(x)\ < 



Cp{xf 



for < a < 2m. 
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We denote by \£ t ne unique function in the homogeneous Sobolev space 
W™ equaling 1 at £, vanishing on and having minimal semi-norm 

I ' lwy i (R d ) over an functions sharing these properties. That is, \£ ls the La- 
grange function for surface spline interpolation centered at £ G S. The La- 
grange functions permit us to rewrite the interpolant as Is/ = X^es /(£)%£• 

Our goal is to bound the operator norm of the interpolation operator on 
weighted spaces. By a a simple application of Holder's inequality we 
observe that 



Infix) 




p(x) s 





p{x) 



< 



/(£) 




Assuming a 'slow growth assumption' introduced in [8] and defined in the fol- 



lowing section, it suffices to bound L a%p := sup^gn X^ges \^ + -pj^y 

for any nonnegative a = s(l — e). For a target function / having support 
compactly contained in f2 (as we assume throughout this article) the sum- 
mation is over the subset S/, so we need only bound 



(2.2) 



We note that, when cr = (or when p is constant), this is precisely the 
usual Lebesgue constant, and for this reason, we call C a)P the penalized 
Lebesgue constant. We approach this by estimating the Lagrange functions 
associated with £ slightly removed from the boundary of Q. The result we 
obtain, under assumptions placed on the centers, detailed in the subsequent 
sections, (Theorem I4.3|) is 



\Xt(*)\<C[l 



\x-(\ 

p(0 



exp 



-A 



\x-Z\ 

p(0 



for all x in the largest ball centered at £ that is contained in f2. 



3 Local Density Functions 

We begin by developing the notion of a local density (LD) parameter. Such 
parameters are key to the results of [4J, although we follow the approach 
taken in [8j. Initially, we can think of the LD as a function on which, at 
a G fi, is the minimal radius needed to capture a ET-stable local polynomial 
reproduction of order £. 
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Definition 3.1 (Local Density Parameter). Given a set of centers 3 C M. d , 
a local density parameter ( LD ) p : W 1 — > R + is a function with an associated 
local polynomial reproduction a : (£, a) i— > a(£, a) : S x R d — )• R of precision 

& 

(Support) For |£ — ar| > /?(a), a(£,a:) = 0. 

(Precision) For all p G ID. we have X^eH a (£' q )p(0 = P( a )- 

(Stability) There is K > such that X^geH l a (£; Q )l < ^ f° r au a - 

This definition is not quite sufficient for our purposes, and we impose an 
extra condition of global compatibility to ensure that the LD does not grow 
too rapidly. This is the approach taken in [8] and which is at the heart of 
DeVore and Ron's error estimates. 

Definition 3.2 (Slow Growth). We say p exhibits 1 — e slow growth if 

p(a)<p(x)(l+ l -^^j , (3.1) 

for every x and a in £1. 

The function p : Q — > R+ exhibits self-majorization of order r if there is 
a constant C sm > so that for every x and y in f2, we have 

P(y) > C smP (x) (l + ^y^) • (3-2) 

In [8] it has been shown that these conditions, (|3.ip and (|3,2p . are actually 
equivalent assumptions, provided r+1 = 1/e. Ifp satisfies 1 — e slow growth, 
then it satisfies self-majorization of order r, with constant C sm depending 
on e only. 

Norming Sets and the Zeros Lemma 

The norming set concept has become an essential tool in scattered data 
approximation since being introduced in this context in [TT]. In our setting, 
it allows us to transfer the norm of a polynomial on a domain (such as 
a ball) to the norm of its restriction to a subset. A set T is a norming set 
for a ball B with norming constant k if for polynomials p S II^ we have 

||p|U„(B) < kIIpItII/^T)- (3-3) 

A set of centers associated with a local density function provides precisely 
this: 



S 



Lemma 3.3. Suppose S has a LD p with precision t, stability K that satis- 
fies slow growth. Then E x := EnB(x, 3p(x)) is a norming set for polynomials 
p 6llf in the ball B(x,3p(x)). 

Proof. The maximum value of p(y) for y £ B(x,p(x)) is bounded by 2p(x), 
so for every y £ B(x,p(x)) there exists a local polynomial reproduction 
a (y,0 °f order i having support inside B(y,2p(x)) C B(x,3p(x)). Thus 



\p(y)\ < 



«6S 



□ 



This brings us to the main result of this section: the zeros lemma, which 
states that the norm of a function in a Sobolev space, vanishing on a suffi- 
ciently large set, is bounded by a Sobolev seminorm of that function. 

Lemma 3.4 (Zeros). Suppose x £ R d , m > d/2,andr > 0. //T C B(x,3r) 
is a norming set for p £ H m -i with norming constant k, then there is a 
constant C depending on m, k and d such that if f £ W™(B(x, 3r)) vanishes 
on T, then 

^ u ^W%{B{x,r)) - u 'w 2 m (B(x,3r)) 

k<m 

Proof We begin by assuming r = 1, x = 0, and writing B\ = B(0, 1) and 
£>3 = B(0,3). For any / £ W2 1 (B^), there is a polynomial p of degree at 
most m — 1 so that ||/ — p||u/ 2 m (s 3 ) < C|/lw 2 m (B 3 ) (this is sometimes called 
the Bramble-Hilbert Lemma). We can estimate ||p||vi/ 2 m (Bi) by its L 00 (Bi) 
norm, with a constant depending only on d and m, and, consequently by 
the £oo norm incurring a constant k, 

||p||w3»(Bi) < C\\p\\ Loo(Bl) < CK\\p\ r \\ tgo (jy 

By our assumption, / vanishes on T, so there is no cost in subtracting its 
restriction from p\ r . Thus, we get 

||p||wy(Bi) < C||(p-/)| t II<oo(T) < C|b-/IUoo(B3) ^ C\\p - f\\ w ™ {B:i) . 

The last inequality is an application of the Sobolev embedding theorem. It 
follows that 

ll/llw 2 m (Bi) < 11/ -pWw^iB-i) + IMIw 2 m (Bi) < (1 + C)l/lvK 2 m (B 3 )- 

The result for general r and x holds by applying an affine change of variable, 
y 1 y x ~\~ ry. □ 
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Remark 3.5. Property (j3.3|) can be modified to treat slightly smaller subsets 
- a fact we use later. In particular, if ()3.3[) holds for ILj, then a similar result 
holds for polynomials of lower degree by removing one point from T. Let 
T' = T \ {£}, and consider p G He-i- We wish to show that 



IpWl^b) < kC £ 



P\ 



Without loss of generality, we may consider B = B(0, 1) and T C B(0,3), 
since the result for general radii follows by rescaling. 

There is a constant q depending only on £ (but independent of £) so 
that \\p\\ Lao (B) < Ci\\p\\ Loc (n) where Q = B \ B(£,, 1/2). By compactness, we 
can choose y G cl(B) so that \y — £| > 1/2 and c^\p{y)\ > \\p\\b- Consider 
P G He, where P(z) := ((z — £), ^jE^)p( z )- Then 

IIpIIloo(b) < 2c e\v - £1 \p(y)\ = 2c e \P(y)\ < 2c £ k \\P\r\\ ioo( r) ■ 

The last inequality follows by applying (|3.3p to P. 

Observe that -P(£) = 0, so H-Plxll^m = \\P It'II^^y'V Since |£ — C\ < 6 
for every ( G T', it follows that \\p\\ Loo ( B ) < ^wll-PH^T') < 12 c e K\\p\\ £oc{r) . 



4 Lagrange function estimates 

To show that the Lagrange function is rapidly decaying, we make use of 
a type of discrete Gronwall lemma. The key estimate is to show that the 
bulk of the tail of an exponentially decaying function can be captured over 
a sufficiently wide annulus. For an increasing function R : M + -4 M+ we 
express this as: 

^ X ^W^(B(^R(t))\B^,R(t-l))) ~ ^ ~ ^I^W 2 m (lR d \B(e,-R(t-l))) ' 

for < /i < 1 independent of t. This implies the complementary inequality, 

^W^\B(t,R{t))) ~ ' Whkh ^ iterat6d lt ~ 1] 

times to obtain 

This is the type of argument made in Lemma l4.1l and Lemma l4.2i 

Lemma 4.1. Assume that Set! has an LD p, with precision m — 1 satis- 
fying 1 — e slow growth. Suppose, furthermore, that dist(^,(9il) < r$. Then 
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there exists a constant p £ (0,1) (depending on e,m,d and the stability 
constant K of p) so that for sufficiently small p(£), i.e., for p(£) < 7 := 
(er - e )^ ; and 3 < t < (r /p(£)) e , 

' X ^'w 2 m (R d \B(c,p(5)t 1 /«)) - ^l X ?lvy 2 m (]R d \B(5, / 3(0[t-3] 1 /£))- 

Proof. Without loss of generality, we assume £ = and write po := p(£) 
and x := X£- With 1 + r = 1/e, construct a C°° cutoff 0, vanishing outside 
B((t — l) 1+r po) and equaling 1 on B((t — 2) 1+r po). That is, <f> is piecewise 
constant outside of the annulus 

A m :=A((t-l) 1+T p ,w(t)) 

which has outer radius (t— l) 1+r po and width w := [(t— 1) 1+T — (i — 2) 1+T ] po- 
We note that (r + l)(t - 2) T p < w(t) < (r + l)(t - l) T p < 1, for p Q less 
than 7. It follows from the chain rule that 



Halloo <C(a)(w)-\ a \. 
By the variational property of x, we have that 



(4.1) 



\x\ws 



W / 2 m (R d ) 



< Ix 



+ 1^x1 



where we express the ball interior to A% n by the compact notation: 

B:=B(p (t-2) 1+T ). 

Our goal for the remainder of the proof is to estimate the quantity \<fix\\ 
This takes the form of two estimates. The first estimate is concerned with 
the effect of multiplying by the cutoff, which can be easily measured using 
Leibniz's rule in conjunction with (|4.ip . This is independent of the radius 
t, and allows us to control the seminorm in terms of a biased Sobolev norm 
of x- The second estimate will control this biased Sobolev norm via parti- 
tioning the annulus into small balls where Lemma 13.41 can be applied. 



Estimate 1 The desired result is to control \(/>x\ 



by a weighted 



combination of Sobolev seminorms of \ on A{ n . We observe first of all that 
we can simplify matters by applying the product rule: 

2 



|0X| 2 



W 2 m (A„) 



E 

|o;|=m 



ml 



ft! 



E 



D^ X (x)D 



dx 



< 



ci E E 

|a|<m /3<a 



D p X {x)D a - f3 (, 



dx. 
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The constant C\ = C\ (m, d) is seen to depend only on m and d. We continue 
with the estimate, now observing that we can remove the factors D a "^(i) 
by flUE]). Thus, 



wxftw-w s c > E E f (») 2(I " H °" 

|a|<m/3<a - 4 " 1 



2 



dx 



< GEW'-'IXI^^,. (4.2) 



The first inequality follows by applying (|4.ip and the observation that any 
two multi-indices f3 < a satisfy \a — f3\ = \a\ — \/3\. The second inequality is 
obtained by rearranging terms. 

Estimate 2 We are nearly in a position to apply Lemma 13.41 Cover Ai n 
with a sequence of balls (Bj)j^j such that 

• each ball is of radius 3pot T , 

• each ball has its center in Ai n and, thus, Bj C A ou t where 

A out := B(po[(t - 1) 1+T + 3t T }) \ B(p [(t - 2) 1+T - 3t r ]), 

• the balls have a finite intersection property: every x £ Aout is in at 
most Nd balls Bj, with Nd depending only on the spatial dimension 
and not on i, po,T, etc. 

Writing the center of each ball Bj as Cj, we observe that the slow growth 
assumption implies that the density p(cj) is not greater than p (l + (t — 

l) 1+r ) < p t T , and LemmaO may be applied on each ball. The expres- 
sion from fj4.2j) can be bounded by corresponding norms carried by the balls 

Bj- 

E(*) 2, " m, i^<A„, s EEW !M wi,;w 

k<m k<.mjaj 

2(k-m) 



<- EE(^ +1 )(i f ) 7 



I X I wfc 



W${B 3 ) 



j&J k<m 



The second inequality follows from the observation that, because t > 3, 
the width of the annulus, which satisfies w > po(r + l)(i — 2) T , is greater 
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than pq(t + 1) (|t) 7 . The third inequality follows from Lemma 13.41 with 

'II 
,3> 



c = (r + 1) (^) T+1 and k the constant from Lemma [3.41 The final inequality 



follows from the finite intersection property, and C = NdK 2 c 2m . Combining 
this with ( 14.21 ). we see that 

|2 , 1 2 ^ I j.,,|2 ^ u,|2 



with C 4 = C 4 (m,d) = C 3 C. 

It follows from their definitions that B = B (p (t - 2) 1+T ) C B (p t 1+T ) 
which implies B^ (po^ 1+r ) CS''. In turn, we can rewrite the annulus as the 
difference of two punctured planes A ou t C (po(t — 3) 1+T ) \B (/?oi 1+T ) ■ 
Applying this to fj4.3j) we get 

Mw 2 " J (B c (p t 1+T )) - Ci i} X \w,^(B c (po(t-3) 1 +^)) ~ lxlwf 1 (B c (p * 1+T ))) 
and, consequently: 

. . . vm , _ 

with " : = 7& < L D 

Iterating this, we show that the seminorm of the characteristic function 
on a punctured space decays exponentially. 



Lemma 4.2. Under the assumptions of Lemma \4-l\ there are constants 
C, v > depending only on e, d, m and K so that 



Ixdw2 m (-B C (£,T)) - ^\Xi\w^{K. d ) ex P 



T 



AO, 

holds for < 7 (the constant from Lemma \4. 1\ ) and 3 T+1 /?(^) < T < vq. 

Proof. As before we take £ = and write /?(£) = p(0) = pQ. Let T = pQt T+1 , 
and let N be the greatest integer less than | (pB') ~ 3 ■ This represents 
the number of times Lemma 14.11 can be iterated. Thus, we see that 

\x\w^(B^{T)) = \X\w^(B c (p f r + 1 )) 

< l 1 \x\w^(B^(po(t-3y+ 1 )) 



The result follows with e v = pa and C = p l . □ 
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Let := dist(£, H \ {£}) be the 'separation distance' at £. By com- 
paring the seminorm of the Lagrange function to that of a suitably scaled 
bump, it is possible to estimate \xi\m in terms of <?(£), the distance from £ 
to its nearest neighbor. On the other hand, by applying the Sobolev embed- 
ding lemma for balls, in conjunction with the zeros lemma, Lemma |3.4( it is 
possible to obtain pointwise estimates for Xi- This is done in the following 
theorem. 

Theorem 4.3. Let HcU denote a finite set of centers with associated LD 
p of precision I > m and stability K, satisfying 1 — e slow growth. Then 
there is A > depending on e,m,d and K (specifically, A := v/2 e , with v 
from Lemma \4-*fy so that for every £ G S, with dist(£, dQ) < ro 



i«Mi<c(fir i/2 (i + ^f) «p 



> /min(|x - £|,r ) 



with s = (1 — e)(m — d/2). 

Proof. We may compare the Lagrange function X£ to another (worse) inter- 
polant a (jj^j)i to make the initial observation that 

\xt\m < (<?(£)) "VL- 

The Sobolev embedding theorem implies that |u(0)| < C||«[|^-m(B(o,l))- Set- 
ting u(y) = X( { x + p( x )y) > we have 

i»wi<cIE"W 2 "i«i^ (b(i , (>(i))) v 

V k<m 



By Remark 13.51 (S n -B(x,3/?(x)) \ {£} forms a norming set for LT m in 
B(x,p(x)). Thus, Lemma [33] implies 



1/2 



|xe(x)| < cp(xr- d / 2 ^ P (x) 2 ( fc - m )|xd 

\ k<m 



i fe<m 

\ m— d/2 

.9(0' J 

l-e 

result for /?(£) > 7, since in this case the factor exp 
C, a constant depending only on ro,e and A. 



Applying the slow growth assumption p(x) < (l + gives the 

\ / min(|a-g|,r ) , 

A| — m — ' ' 
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This also proves the theorem for the case p(£) < 7 and |£ — x\ < k e p{^), 
where 

k e := e 1 /^ 1 -^ 6 . 

Suppose p(£) < 7 and |£ — x\ > fc e p(£). For T = dist(£, x) in the range 
3 T+1 p(£) < T < ro, we again consider centers in the ball, -B(x, 3p(x)), where 
Lemma 13.41 can be applied. By a similar argument to the one employed 
already, we have 



< C{p{x)/q{i)) 



m-d/2 



cxp 



))) 

T - 3p{x) 
P(0 



Writing T := ^/J, , the bracketed statement in the exponent can be written 



as 



< [- I /(r-3(2r) 1 - e ) ei 



where we invoke slow growth and the fact that T > 1 to estimate < 
(2r) 1 ~ £ . Observe that T > k e = 6 x / e ■ 2^^, implies that T - 3(2r) 1 " e > 
r/2, so [-v (v - 3j^)] < [-A(r) e ] with A = v/2 e , and the theorem fol- 
lows. □ 



5 Penalized Lebesgue Constant Estimates and Main 
Results 

To apply the result of the previous section, we need to estimate the con- 
tribution at x from every Lagrange function X£- Despite their fast decay, 
we will never succeed if we do not control the number of centers jp£. Nor- 
mally this would be accomplished with a constraint on the spacing of H, a 
so-called quasi-uniformity condition, which would permit control of #H by 
a function of the global density h. In our setting this is out of the ques- 
tion. Rather, we make the (mostly harmless) assumption that the distance 
between two nearby centers is controlled from below by the local density: 
g(£) = min^gn^^ |£ — CI > We call such an assumption weak quasi- 

uniformity. The ratio of weak quasi-uniformity is cq = max^= p(£)/<?(0- 

Theorem 5.1. Assume f E C C (Q) satisfies dist(supp(/), diij =: vq. Fur- 
thermore, assume 5 C 0, has local density function p of precision £ > m 
and stability K , satisfying 1 — e slow growth, and weak quasi-uniformity, 
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with ratio cq. Then the penalized Lebesgue constant C p ^ a is bounded, with a 
constant C(tt,cr,d,e,co, K,ro) depending only on £l,o~, e,co,K and ro (and 
not on p). 

Proof. Set o"i = cr+(l— e)(m— d/2). Theorem 14. 31 permits us to estimate the 
penalized Lebesgue constant as £ PjCT < C Yl^eS + / ex P 
I + II, where 

' - § (-^r)"-h(^))l 

|a'-g|>ro 
|a:-£|<ro 

To treat /, we decompose the sum into subsets of S based on the size of 

P- 

E k := {£ G S | < p(0 < 2- k }. 

Since f2 is bounded, p is bounded above by r\ and, hence, there is k$ G Z so 
that p(£) < 2- fe °, for all f. By weak quasi-uniformity, there is k so that for 
any k G Z, 2-( fc+K ) < g(f) for all and < C d |ft|2- d ( fc+K ), which 

implies 

oo 

1 ^ EE (l + n^'exp 
k=k £e= fc 
oo 

< C \n\2( k+K)d 2- kai exp 
fe=fco 

with < b := Ar^. Because e > 0, the final sum is bounded, with a sum 
depending only on ro, e, r\, Co, |fJ|, eZ and a. 

To treat II, we employ subsets of E defined according to the distance 
from x 

X k := U G S | 2 k p(x) < \x - £| < 2 fc+1 p(x)}. 

Slow growth and self-majorization allow us to estimate p(£) from above 
and below on Invoking the self-majorization of p, we see that p(£) > 
p(x)2- T ( fc+2 )-- ?0 for £ e X k , and for some number jo depending only on 
e for which 2~ J0 < C sm . This has two consequences. First, Kj^ < 



min |as— £|,ro 



-A (r 2 fc ) 
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2 jo+r(fc+2)l£_5i < 2 i 1+k ( T+1 \ which implies the estimate 

pkO ) 

Second, quasi- uniformity ensures that <?(£) > cop(x)2- T( - k+ V-i° , which al- 
lows us to estimate #X k by {2 k+1 p{x)) d /(c p(x)2- T( - k+ ^-^) d . Hence < 
^2 rffc ( r + 1 ) with a constant depending on e, (i and Co- The greatest value of p in 
X k is not larger than p(x)2 ( - k+2 ^ 1 - t \ It follows that J|^l > frl^gE^ -1 ) > 
b2 ke , for b > depending on e. Thus 



exp 



-A 



\x-£\ 

p(0 



< C^2*( T+1 )2 /c ( T+1 )' J1 exp 



-y2- fee2 



fc=o 



with A' = A6 e > 0. It follows that II is bounded by a constant depending 
only on e, d, and <r. □ 

The following corollary applies Lebesgue's lemma to the local approxi- 
mation results of [4J. Of special interest is the fact that the schemes used by 
DeVore and Ron to produce local results are quite abstract and do not lend 
themselves to direct implementation. For functions of full smoothness 2m, 
the approximant is determined by modifying an integral representation in- 
volving the surface spline. This is worsened when functions of less than full 
smoothness are treated. In this case a complicated K- functional argument 
is invoked in order to get the result. Interpolation, in contrast, is directly 
implementable and universal in the sense that it requires only the data to 
work, and makes no extra requirements of the target function. 

Corollary 5.2. Assume that H is a set of centers in having a local den- 
sity function p of precision £ greater than 2m + 1 — d, with 1 — e slow 
growth, weak quasi-uniformity and, furthermore, that f £ C(f2) satisfies 
dist(supp(/),afi) > r . If f € A s , with A s = C 2m {9.) when s = 2m or 
^00 00 0^) w hen < s < 2m, then 

\f(x)-I 3 f(x)\<Cp(x) s ||/|U. (n) 

with a constant depending on f2, m, e, cq and tq. 
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Proof. We must show that there is a constant C so that, for any target 
function / 6 A s and LD p satisfying the conditions of the corollary, there 
exists Sf- £ S- so that 



< C. 



(5.1) 



In [3], an operator To : C,?" 1 (M ) — > span^ gS <j)(- — £) is constructed mapping 
compactly supported functions to linear combinations of the RBF 



Tsf 



A m f(a)a(£,a)da 



0, 



with the desired approximation property: i.e., sjo = Tsf satisfies (|5.ip . 
For functions of lower smoothness (i.e., / S 00 ), (|5.ip holds by first ap- 
proximating / by a smooth function g satisfying the properties |A m (7(x)| < 
C7/j(x)*- 2m ||/||B» 0i00 and \f{x)-g{x)\ < C7^(z)*||/||b», i00 and then selecting 
s/n '■= T^g. In either case, the approximant s/n has coefficients determined 
by an expression of the form = J Rd A m F(a)a(£, a)da, where F has com- 
pact support. Utilizing the polynomial reproduction of a, we see that, for 

P ^ n min (£ 2m-l)i 

V / A m F{a)a{£ : a)p{£)da= I A m F(a)p{a)da = 0, 



by applying Green's identity 



A m F(a) p(a)da 



'B(0,R) 
2 m- 1 



F(a)A m p(a)da 



A j+ iF(a)A 2m -jp(a)d<T(a) 



dB(0,R) 



for a ball B(0,R) of appropriate radius (the boundary operators Aj are 
Aj := A^ -1 )/ 2 or Aj := D n A^^ 2 ^ 2 for odd and even j, respectively). Thus, 
sj = 6 and the corollary follows from Lebesgue's lemma. □ 
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